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O ! Abstract 

o 

\ We consider the effects of photon diffusion on the small-angle microwave background anisotro- 

pics due to active source models. We find that fluctuations created just before the time of 
last scattering allow anisotropy to be created on scales much smaller than allowed by standard 
Silk damping. Using simple models for string and texture structure functions as examples, we 
illustrate differences in the angular power spectrum of the intrinsic and Doppler components of 
Q\ ■ the anisotropy on scales of order a few arcminutes. In particular, we find that the Doppler peak 
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heights are modified by 10-50% and small-angle fall-off is power law rather than exponential. 



1. Introduction 

Accurate measurements of the anisotropy in the cosmic microwave background radiation are a 
powerful probe of modern cosmological models since they represent a snapshot of the universe 
just before it became transparent to electromagnetic radiation about 400,000 years after the 
Big Bang. The standard picture is that the initial seed fluctuations were created close to 
the Planck epoch by quantum effects, during a period of super-luminal expansion known as 
inflation. These perturbations have been termed as 'passive' [1] since they were created during 
an early epoch and follow a linear deterministic evolution until last scattering, where they leave 
an imprint in the microwave background. 



Email : r.battye@ic.ac.uk 

Paper submitted to Physical Review D. 



1 



There are, however, a potentially much larger group of theories, known as 'active' source 
models [2,3], which could lead to the creation of microwave anisotropies and the formation 
of large-scale structure. The best motivated of these are those associated with topological 
defects [4,5,6] formed during cosmological phase transitions, such as cosmic strings [7,8] and 
textures [9]; their generic feature being that they create perturbations over a very wide range 
of scales at all times, from their formation to the present day. 

The large-angle (^1°) anisotropies due to defects have been shown to be compatible with 
those detected by COBE since they lead to a near scale-invariant spectrum of perturbations on 
the relevant scales (see, for example, [10]), while little quantitative work exists on the formation 
of large-scale structure due to calculational difficulties in modelling non-linear effects. The 
small-angle (^$1°) microwave background anisotropies provide an area where the predictions 
from passive and active theories can be tested since they are generically very different. A 
number of works [1,2,3,11,12] have predicted shifts in the positions of the so called Doppler or 
Sakharov peaks which occur largely due to the differences in the tight coupling solutions. Here, 
we shall discuss the effects of photon diffusion or Silk damping [13] which have until now been 
ignored or introduced ad hoc. 

The crucial period for understanding these effects is that just before last scattering. During 
this epoch the acoustic oscillations in the photon-baryon fluid are damped by the increasing 
mean free path of the photons. It is not difficult to see that fluctuations created after the onset 
of this regime will receive less damping than those created before it. In particular, those created 
just before the time of last scattering will receive virtually no damping at all. If perturbations 
are created on all scales above the defect size, as is thought to be the case for defect models, 
then it will be possible for anisotropy on small angular scales to remain to the present day. 
We will see that it is not sufficient to model these effects with a simple multiplication by an 
exponential suppression factor across all scales. Rather it requires careful consideration of the 
time at which fluctuations are created. Simple estimates will show that there are potentially 
important modifications to peak heights and also power-law suppression at small angular scales, 
rather than exponential. 

2 Analytic formalism 

We shall use the analytic formalism developed in refs. [14,15] to describe the small-angle mi- 
crowave background anisotropies for passive perturbations. There are, however, shortcomings 
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in the way Silk damping is modelled which will lead to discrepancies at small angular scales 
when applied to active sources. The starting point for this treatment is the first-order collisional 
Boltzmann equation in Fourier space ignoring the effects of polarisation, 

e + zM© + *) = + « (Go - e - ±p 2 (fj)e 2 - ip^v*) . (1) 

In this equation, Q(k, 77, ]i) = J2(—i) l Qi(k, i])Pi((i) is the Newtonian temperature perturbation, 
77 is conformal time, k is the wavenumber, Pi(fi) is the Legendre polynomial associated with 
angular variable ]i = cos 6, Vb is the baryon velocity, k is the differential optical depth due 
to Thomson scattering and ^ , $ are the gauge invariant potentials [16] assumed to be char- 
acterised by an external source. One can solve this equation for © ; (/ > 1), if one allows the 
intrinsic (or monopole) temperature anisotropy and the baryon velocity to act as extra source 
terms. The solution at the present day is then [17] 

Vo , s rVo 

dr] e ~ K ^^ o) - $J ji [kAr]] + J dr] ke~ K ^ m) - $) j t [kAr]] 
+ J d V ke- K ^vo) (e ojl [( k Ar,] + ^—^Vbj^kAr]] - -J—V h j l+1 [kAr,] 

where Go = ©0 + Ar] = r]o — r], ^(771, 772) = dr]k(r]) and ji(x) is a spherical Bessel 
function. This expression is slightly different to previous formulations, in that there are the 
separate terms 60 + $ and * — <E», as opposed to just ©0 + The reason for this is that, as we 
will see, it is only O + <&, henceforth described as the intrinsic anisotropy, which is modified 
by Silk damping. 

The function known as the 'visibility function' and is sharply peaked around 

the time of last scattering 77*, while e - ^ 77 ' 770 ) is zero for 77 << 77* and effectively constant 
77 >> 77*. Therefore, one can see that the anisotropy is essentially due to three effects. The 
first sums up all contributions due to the motion of the sources after the time of last scattering 
and the second is just the difference of the two potentials at last scattering. The combination 
of these two terms, known as the Sachs- Wolfe effect [18], dominates the anisotropy on large 
angular scales and possibly again at very small (^$1') angular scales where, for example, the 
Kaiser-Stebbins effect [19] is prevalent for cosmic strings. The final contribution is due to 
the damped acoustic oscillations of the intrinsic anisotropy and the Doppler effect around last 
scattering and it is this which we shall consider in this paper. This is not to say that the other 

3 



®i{k,r]o) 
21 + 1 




contributions will have not effect over the range of angular scales considered, rather that there 
are other effects — likely to be of similar amplitude — which need to be taken into account. 

It is clear that in order to calculate the small-angle anisotropy, we must first calculate 60 
and Vb around the time of last scattering and then deduce 6; using (2). To do this we decouple 
the Boltzmann equation into its multipole moments, 

O = -* - ifc6i , 9! = fc(9 + *) - |«9 2 - A(6i - v b ) , 

for / > 2. As it stands (3) is not a complete system of equations, which can be rectified by 
introducing the Euler equation for the baryon velocity 

v h = --v h + kv + ^-(e 1 -v h ), (4) 

a R 

where R = 3pt>/4p 7 is the baryon to photon ratio normalised to be 3/4 at photon-baryon 
equality. Thomson scattering is highly efficient at early times with k large, therefore expanding 
(3) by ignoring O(kfk) yields the equations for tight-coupling regime. The effects of Silk 
damping can be investigated by expanding to higher order. If one ignores 0(k 2 /k 2 ) and also 
0(R 2 ), the equation for the intrinsic anisotropy is 

* + (its + ^tt*) g » + - - 1* 2 (its - *) (5) 

and the baryon velocity is given by the dipole temperature anisotropy, V\> = ©i = — 3©o//c. 
One can solve (5) for large k using the WKB approximation and Green's method, 

[1 + R( V )} 1/4 §0(77) = e - fe2 / fc ^>°) © O (0)cosfcr s (r ? ) + ^ (§ (0) + ±£(0)e (0)) smkr s (rj) 

+ ^ f V dq' [1 + R(v')f /4 e - k2 / k * ( ™') S in [kr s (r]) - kr s (rj')} H{rf) , 
k Jo 

(6) 

where r s {rj) is the sound horizon distance and /c^~ 1 (?72, Vi) is the Silk damping length, below 
which photon diffusion removes anisotropy, 

1 p drf 2 4 /"» drf 1 m 

V3 Jo V 1 + #(*7 ) 27 ) 1 + ) 



This solution has two parts: the transient whose amplitude and phase is dependent on the 
initial conditions and the particular integral due to the forcing term. At the simplest level, the 
transient solution corresponds to the contribution from passive perturbations, while calculating 
the particular integral represents the effects of active perturbations by a convolution of the 
source term and the oscillatory terms which represent the acoustic waves with the appropriate 
damping coefficient"!". The effect of Silk damping is clearly different for the passive and active 
perturbations contrary to the approach of refs. [14,15], since the damping length for active 
perturbations depends on when the perturbation is created. 

In order to investigate the consequences of this effect on the microwave background, we 
shall make a number of simplifying assumptions. These assumptions allow simple calculations 
to be made, but the underlying physical principles which we want to illustrate will remain. First, 
we assume R ~ 0, which will be a good approximation in the low baryon content universes 
predicted by Big-Bang nucleosynthesis, since R(r]*) ~ 30fit>/i 2 << 1 t- Also we assume only 
active perturbations, by ignoring the transient part of the solution. Under these assumptions, 
the sound speed of the photon-baryon fluid is 1/ a/3 and 
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where H( V ) = k 2 ($ - *)/3 and K^m) = (4/27) drf'/k{rf'). 

At this stage, we must attempt to specify the nature of the source. This can be done 
using a structure function [21], a two-point correlation function of the source potentials, which 
determines the wavenumbers that are perturbed at a particular time. A simple example is 
x 4 F(x) ~ S(x — x c ), for which a single wavenumber k = x c jr\ will be perturbed at conformal 



t It should be noted that passive perturbations also create perturbations at horizon crossing and hence the particular 
integral will also contribute. However, it can be seen that a source creating perturbations on lengthscales close to the 
horizon size will lead to a contribution close to being in phase with the initial conditions. 

| This approximation ignores the effects of baryon drag [20] , which lead to a power- law fall-off in the power spectrum 
at small angular scales. Since this effect is also due to the non-trivial time dependence of the gravitational potential 
it is likely to be of a similar size to the power-law fall-off discussed here, although the exact relative size will depend 
crucially on the baryon density. 
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Figure 1(a): The effects of the modified damping formalism compared to the tight coupling solution 
and the standard Silk damping formalism for coherent statistics. We plot the power spectrum of intrinsic 
component of the anisotropy on a log- linear scale. On the left, the tight coupling solution for the cosmic 
string (bottom) and texture structure (top) functions, in the centre the standard damping formalism and 
on the right the modified formalism. It is clear that the modified formalism allows much more power to 
remain on very small scales, with a power law fall-off. 

time 77. More realistic examples of active perturbations will generate perturbations over a range 
of wavelengths. Of particular interest in this context is the asymptotic fall-off of the structure 
function: for small x causality [22] suggests that F(x) ~ const, whereas for large x the fall-off 
is given by the characteristic shape of the anisotropies created by a particular model. Cosmic 
strings form line-like patterns [19] leading to x 4 F(x) ~ 0(x~ 2 ), whereas textures form spots 
implying that x 4 F(x) ~ 0(x~ 4 ). For definiteness, we use F(x) = (x 2 + a 2 )~ 2 (l — bx + cx 2 )~ n 
with a = 2.45, b = 0.12, c = 0.006 and n=l or 2. If n = 1 the structure function may model a 
cosmic string network [3], while n = 2 will be described as a texture model, although in reality 
the precise values of a, b and c are likely to be somewhat different. These models have been 
chosen to exemplify the possible effects of the modified formalism, although they do have some 
relevance for the realistic models which one eventually wants to describe. 

3 Semi-analytic estimates 

For the purposes of this paper, we first concentrate on calculating the power spectrum (denoted 
by V) of the intrinsic and Doppler contributions to the anisotropies, exhibiting the differences 
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Figure 1(b): The effects of the modified damping formalism (dotted line) compared to the standard Silk 
damping formalism (solid line) for coherent statistics. We plot the intrinsic anisotropy on a linear scale for, 
on the left the cosmic string model and on the right the texture model. This exhibits the differences in the 
peak height for the secondary peaks. Note that the scale starts at X* — 8, which ignores the first peak. 
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Figure 2: The intrinsic component of the anisotropy for incoherent statistics using the same conventions 
as fig. 1(a). Once again the modified damping formalism creates more power on small scales. 
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between the simple exponential suppression — described as the standard Silk damping formal- 
ism — and the more accurate modified damping formalism described in the previous section. 
This allows the study of the separate physical effects without confusion, within a calculationally 
more simple and intuitive framework. 

In order to do this one must specify the type of statistics the source satisfies. A simple 
approximation often used is to assume that the source is totally coherent. In this case P($ — 
vp) = t] 3 F(kr]) [2,3] and hence the power spectrum of the intrinsic contribution to the anisotropy 
is given by 

k s V[Q + ®] = l 

(9) 

as a function of the dimensionless parameter x* = krj^/y/3. However, if the source is totally 
incoherent then various modifications to our understanding must be made [2]. In this case, 

^ + *=\r"** m ^( x '-%)^{-^(*'^)) ■ <io) 

The difference between these two expressions is that the coherent approximation allows con- 
tributions to the anisotropy created at times very much before the last scattering to cancel 
out, whereas the incoherent approximation prevents any cancellations whatsoever, leading to 
the absence of secondary Doppler peaks [2,3]. Neither of the approximations is likely to be 
a complete representation of the true physical situation on all scales, although one may be 
dominant. Here, we shall present results for both approximations, but the reader should note 
that the texture models are likely to be more coherent than string models. 

One can investigate the effects of damping by evaluating intrinsic anisotropy power spectra 
(9) or (10) for particular structure functions using a numerical integration routine and the simple 
model for the standard reionization history (fi = 1,0b = 0.05, h = 0.5) presented in ref. [14]. 
Figs. l(a),(b) show the results for the coherent approximation. Both models have their first 
peak at x* ~ 4 since that is only dependent on the coefficients a, b and c which have been 
kept constant. If damping is included, differences between the standard and modified spectra 
become noticeable for x* > 10 and dominate for x* ^50; the effect being very much stronger for 
strings since they create more power on the very smallest scales. As fig. 1(b) shows the peak 
heights are modified by an increasing amount as increases with, for example, in the cosmic 
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string model the second, fourth and six peak heights being increased by 13%, 25% and 44% 
respectively. 

The decay of the spectrum is approximately exponential for x* < 50 as in the passive 
case, but with a slightly different coefficient which leads to the modifications in the peak 
heights. However, for x* >50 the behaviour of the spectrum is very different; power law decay 
with virtually no peaks. This distinction can be understood by realising that for x* ^50, a 
substantial proportion of the power is created before the last scattering epoch and hence the 
spectrum experiences almost the full effect of the damping, while for > 50 the power is created 
during the last few oscillations of the radiation fluid, preventing cancellations and hence peaks, 
with the power law decay reflecting the large x fall-off of the structure function. The same 
features are seen in the power spectrum for the texture model, but the modifications to the 
peak height are smaller and the transition to the power law regime is shifted to smaller scales 
(x* ~ 60) since less power is created on the smallest scales. 

Fig. 2 shows the equivalent spectra for the incoherent approximation. The tight coupling 
solutions are now totally devoid of any oscillations since no cancellations have taken place. This 
seems to lead to a further shift in the peak to about x* ~ 8, due to perturbations created at 
earlier times feeding into smaller scales. However, the general picture of the modified formalism 
differing from the standard picture for x* > 10 remains. 

The power spectrum of the Doppler (or dipole) component was also investigated by re- 
placing the sine in (9) and (10) with a cosine. We have already commented that the anisotropy 
on small angular scales is dominated by the last few oscillation of the radiation fluid. Formally 
this corresponds to around the region of the top limit of the integral in (9) and (10). For the 
intrinsic component the integrand is zero at the upper limit, but the situation is very much 
different for the Doppler component which is totally out of phase with that of the intrinsic 
component. It will never be zero in the vicinity of the upper limit of the integral and hence 
the effect is enhanced. Fig. 3 illustrates this for the coherent and incoherent string models 
Of course the Doppler component is suppressed with respect to the intrinsic anisotropy for 
R > 0, but this will have significant implications for the spectrum of polarisation and density 
perturbations produced since they are created by the dipole anisotropy. 

4. Discussion and conclusions 

Until now we have concentrated on the power spectrum of the intrinsic and Doppler anisotropies. 
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Figure 3: The effects of the new damping formalism on the Doppler contribution to the anisotropy. On 
the left, the tight coupling solution for a coherent (bottom) and incoherent (top) cosmic string structure 
function, in the centre the standard damping formalism and on the right the modified formalism. Notice 
that the difference between the standard and modified formalisms is more obvious than for the intrinsic 
contribution. 

In order to calculate the power spectrum of the total anisotropy, we first assume that it is 
dominated by the intrinsic component, which is likely to be a good assumption for angular 
scales between 1° and V — although note the discussion of the increased contribution from the 
dipole above. To deduce an expression similar to those provided in ref. [14] , we must be careful 
to include the effects of cancellation on small scales due to finite width Ary* of the last scattering 
surface. This is represented by the rapidly oscillatory nature of the spherical Bessel function at 
large k. Averaging this effect gives (ji[k(rj - rj)}) m (kAi]*)~ 1 / 2 ji[k(r]o - //*)] for k > (A??*) -1 
and hence 
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[k(Vo-V*)} / 
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drfHijf) sin 
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(V - rf) 



(11) 



where T>(k, rj', 770) is the equivalent of the acoustic visibility function in ref. [14], given by 



rvo 

V(k,rj',rj ) = / dr]ke- K 

Jn' 



(12) 



for k < (A??*) -1 and multiplied by a factor of (/cAry*) -1 / 2 for k > (A??*) -1 . 

Fig. 4 shows a comparison between the standard damping mechanism and the modified 
damping mechanism for a finite width last scattering surface (77* /A?]* ss 6) including the effects 
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Figure 4: The angular power spectrum of the anisotropy due to the (dominant) intrinsic component. On 
the bottom, a sample incoherent cosmic string spectrum for the standard Silk damping approach (dotted 
line) and the modified damping formalism with a finite last scattering surface (unbroken line), which 
includes the effects of cancellation damping. The top illustrates the same spectra for the coherent string 
model. Notice the changes in the peak heights even on this log-linear scale, and the power-law rather than 
exponential tail. 

of cancellation damping for the coherent and incoherent string models. We have assumed that 
the anisotropy is dominated by the intrinsic component and have used the simple relation 
/ = krjo to calculate (AT/T) 2 oc / dkk 2 \Qi(r]o,k)/(2l + 1)| 2 . Although this calculation is 
approximate, it shows the important differences of the power spectra. It is clear that the level 
of anisotropy possible at very small angular scales is much larger with the modified damping 
formalism. 

To summarise, we have presented a simple, intuitive formalism to treat the effects of 
photon diffusion on the CMBR due to active sources. Simple models for cosmic string and 
texture scenarios illustrate that the effects of Silk damping mechanism are somewhat different 
to previously thought for scales smaller than the first acoustic peak, particularly in the case of 
strings. At this stage, we do not claim quantitative accuracy, rather that we have illustrated 
an effect which must be taken into account when making comparison between state of the art 
Boltzmann codes and semi-analytic methods similar to those described here, and that it should 
be included in future work. 



11 



Acknowledgements 

Thanks to Andy Albrecht for impetus and encouragement. I have also benefited from conver- 
sations with Diego Harari, Joao Magueijo, Pedro Ferreira, Rob Caldwell, Neil Turok and Paul 
Shellard. I am supported by PPARC postdoctoral fellowship grant GR/K94799. 

References 

1. Albrecht A., Coulson D., Ferreira P. k Magueijo J. [1996], Phys. Rev. Lett. 76, 1413. 

2. Magueijo J., Albrecht A., Coulson D. k Ferreira P. [1996], Phys. Rev. Lett. 76, 2617. 

3. Magueijo J., Albrecht A., Coulson D. k Ferreira P. [1996], |astro-ph/9605047| , Phys. Rev. D in 



press. 

4. Kibble T.W.B. [1976], J. Phys. A9, 1387. 

5. Vilenkin A. & Shellard E.P.S. [1994], Cosmic strings and other topological defects (Cambridge Uni- 
versity Press). 

6. Hindmarsh M.B. k Kibble T.W.B. [1995], Rep. Prog. Phys. 58, 477. 

7. Zel'dovich Ya.B. [1980], M. N. R. A. S. 192, 663. 

8. Vilenkin A. [1981], Phys. Rev. Lett. 46, 1169. Erratum: Phys. Rev. Lett. 46, 1496. 

9. Tuork N. k Spergel D.N. [1989], Phys. Rev. Lett. 64, 2736. 

10. Allen B., Caldwell R., Shellard E.P.S. , Stebbins A. k Veeraraghavan S. [1996], Phys. Rev. Lett, in 
press. 

11. Durrer R., Gangui A. k Sakellariadou M. [1996], Phys. Rev. Lett. 76, 579. 

12. Crittenden R.G. k Turok N. [1995], Phys. Rev. Lett. 75, 2642. 

13. Silk J. [1968], Ap. J. 151, 459. 

14. Hu W. k Sugiyama N. [1995], Ap. J. 444, 489. 

15. Hu W. k Sugiyama N. [1995], Phys. Rev. D51, 2599. 

16. Bardeen J.M [1980], Phys. Rev. D22, 1882. 

17. Hu W. k Sugiyama, N. [1995], Phys. Rev. D50, 627. 

18. Sachs R.K. k Wolfe A.M [1967], Ap. J. 147, 73. 

19. Kaiser N. k Stebbins A. [1984], Nature 310, 391. 

20. Hu W. k Sugiyama, N. [1995], astro-ph/9602019, Ap. J. in press. 

21. Albrecht A. k Stebbins A. [1992], Phys. Rev. Lett. 68, 2121; ibid 69, 2615. 

22. Robinson J. k Wandelt B. [1996], Phys. Rev. D53, 618. 



12 



